Published in ApJ Letters 

Preprint typeset using MTgX style emulateapj v. 12/14/05 



O 

(N 

< 



o 
u 

o 



> 

m 
O 



X 



CONFIRMATION OF ENHANCED DWARF-SENSITIVE ABSORPTION FEATURES IN THE SPECTRA OF MASSIVE 
ELLIPTICAL GALAXIES: FURTHER EVIDENCE FOR A NON-UNIVERSAL INITIAL MASS FUNCTION 

PlETER G. VAN DOKKUm' AND CHARLIE CONROY" 
Published in ApJ Letters 

ABSTRACT 

We recently found that massive cluster elliptical galaxies have strong Nal A8183,8195 and FeH A9916 Wing- 
Ford band absorption, indicating the presence of a very large population of stars with masses < 0.3 Mq. Here 
we test this result by comparing the elliptical galaxy spectra to those of luminous globular clusters associated 
with M3 1 . These globular clusters have similar metallicities, abundance ratios and ages as massive elliptical 
galaxies but their low dynamical mass-to-light ratios rule out steep stellar initial mass functions (IMFs). From 
high quality Keck spectra we find that the dwarf-sensitive absorption lines in globular clusters are significantly 
weaker than in elliptical galaxies, and consistent with normal IMFs. The differences in the Na I and Wing-Ford 
indices are 0.027 ± 0.007 mag and 0.017 ± 0.006 mag respectively. We directly compare the two classes of 
objects by subtracting the averaged globular cluster spectrum from the averaged elliptical galaxy spectrum. The 
difference spectrum is well fit by the difference between a stellar population synthesis model with a bottom- 
heavy IMF and one with a bottom-light IMF. We speculate that the slope of the IMF may vary with velocity 
dispersion, although it is not yet clear what physical mechanism would be responsible for such a relation. 

Subject headings: cosmology: observations — galaxies: abundances — galaxies: evolution - 
content — stars: mass function 



-galaxies: stellar 



L INTRODUCTION 

The form of the stellar initial mass function (IMF) is one 
of the main uncertainties in the interpretation of observations 
of distant galaxies and in our understanding of the conver- 
sion of gas to stars over cosmic time. Most of the stellar 
mass density of the Universe is in the form of low mass stars 
with masses <C 1 Mq. Despite the large contribution of these 
stars to the total stellar mass of galaxies their contribution to 
the integrated light is small, which makes it difficult to mea- 
sure the form of the IMF in galaxies other than the Milky 
Way. In practise it is usually assumed that the IMF in other 
galaxies and at earlier epochs is the same as in the disk of 
the Milky Way: a power-law with logarithmic slope x = -2.3 
at masses > 1 MQwith a turnover at lower masses (Kroupa 
2001; Chabrier 2003; Bastian, Covey, & Meyer 2010). 

As has long been recognized it is possible to measure the 
contribution of dwarf stars to the integrated light of old stel- 
lar populations from absorption features that have a strong 
gravity dependence (e.g., Spinrad 1962; Cohen 1978; Faber 
& French 1980). There are two absorption features in the red 
that are strong in cool dwarfs but very weak or absent in cool 
giants: the Nal A8183,8195 A doublet (e.g., Schiavon et al. 
1997) and the Wing-Ford FeH molecular band at A9916A 
(e.g.. Wing & Ford 1969; Schiavon, Barbuy, & Bruzual A. 
2000). These features reach depths of 30-40 % and 40-50 % 
respectively in the spectra of individual low mass stars. In 
the integrated light of old stellar populations their predicted 
strengths are of order 1-5 % depending on the form of the 
IMF (see van Dokkum & Conroy 2010, herafter vDC). 

In vDC we presented high quality Keck spectra for eight 
massive galaxies in the Coma and Virgo clusters. Both the 
Na I doublet and the Wing-Ford band are quite strong in these 
spectra, in agreement with previous studies (e.g.. Carter, Vis- 
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vanathan, & Pickles 1986; Hardy & Couture 1988; Couture 
& Hardy 1993). As shown in vDC, the absorption is signif- 
icantly stronger (by a factor of ~ 3) than might be expected 
from a lOGyr old Solar metallicity stellar population with a 
Milky-Way IMF. The best-fitting IMF is more dwarf-enriched 
than even the Salpeter (1955) form, requiring a slope x ^ -3 
down to 0. 1 MQ(compared to x = -2.35 for a Salpeter IMF). 

The main uncertainty in vDC is that we rely on stellar pop- 
ulation synthesis models. The models are based on high- 
quality empirical spectra of individual Solar metallicity stars 
from the Rayner, Cushing, & Vacca (2009) IRTF library of 
cool stars, and use up-to-date isochrones (see Conroy, Gunn, 
& White 2009, vDC, and C. Conroy et al., in preparation). 
However, it is known that the stars in elliptical galaxies are 
enhanced in a-elements compared to stars in the Milky Way 
(Worthey, Faber, & Gonzalez 1992) and this cannot be incor- 
porated explicitly in our modeling. 

In this Letter we test whether the enhancement of dwarf- 
sensitive spectral features in elliptical galaxies persists when 
we compare their spectra to those of globular clusters rather 
than to models. The Andromeda galaxy (M31) has a popu- 
lation of old, metal-rich globular that have similar abundance 
patterns as massive elliptical galaxies (Caldwell et al. 2011). 
However, their dynamical mass-to-light (M/L) ratios are so 
low that they must have normal mass functions - or possibly 
even mass functions that are dwarf-deficient ("bottom-light") 
compared to the Milky Way IMF (Strader et al. 201 1). There- 
fore, if the interpretation in vDC is correct, we should see 
much weaker Na I and FeH absorption than in massive ellip- 
tical galaxies, despite the fact that their ages and abundances 
are very similar. 

2. SAMPLE SELECTION AND OBSERVATIONS 

We observed four bright globular clusters associated with 
M31 for which ages and metallicities had been determined 
in other studies: B143, B147, B163, and B193. The clus- 
ters were selected to be bright and old and to have approxi- 
mately Solar iron abundances (Caldwell et al. 2011). For a 
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Fig. 1. — The five globular clusters that were observed with LRIS. 
The main image is an optical image from the Digital Sky Survey 
and the insets are J band images from 2MASS. Four of the clus- 
ters are among the brightest in M31 and have approximately Solar 
metallicity (Caldwell et al. 2011). The relatively faint cluster B189 
was added to the sample to test how sensitive the conclusions are to 
abundance patterns. 



Chabrier (2003) IMF the total number of stars in these four 
clusters is approximately 10^ (of which ^100 should be AGB 
stars), which is sufficient to properly sample the stellar lumi- 
nosity function. We also observed the very metal-rich cluster 
B189 to assess the sensitivity of the FeH and Nal features 
to changes in abundance patterns (see below). The locations 
of the clusters within M31 are indicated in Fig. [T] In Fig. |2] 
their abundances and ages are compared to those of the ellip- 
tical galaxies of vDC. The globular cluster data come from 
Caldwell et al. (2011) and R. Schiavon et al., in prep. Data 
for the four Virgo galaxies were obtained from Trager et al. 
(2000), for an aperture (which is similar to the effective 
aperture of our measurements of Virgo galaxies). The Coma 
measurements are from Thomas et al. (2005) (NGC 4840 and 
NGC4926) and from Harrison et al. (2010) (NGC 4889 and 
IC 3976). All measurements should be on approximately the 
same system, although systematic differences between stud- 
ies cannot be ruled out. The four luminous M31 globular 
clusters have similar ages and abundances as the elliptical 
galaxies. The cluster B189 has a higher iron abundance and 
a-enhancement than most of the ellipticals. Furthermore, it 
is extremely enhanced in C and N, and it has strong NaD ab- 
sorption (R. Schiavon, priv. comm.). 

The globular clusters were observed on December 3, 2010 
with the red arm of the Low Resolution and Imaging Spec- 
trometer (LRIS; Oke et al. 1995) on the Keck I telescope. We 
used the 600 1 mm"' grating blazed at 1 /im. The slit width of 
0"7 gives a spectral resolution of cr = 1.2 A at 9500 A, corre- 
sponding to a velocity resolution of « 40kms"'. Each clus- 
ter was dithered along the slit in a series of four 300 s expo- 
sures. Data reduction followed standard procedures for long- 
slit spectroscopy. Owing to the fully-depleted LBNL CCDs of 
LRIS fringing is not a concern. We corrected for atmospheric 
absorption in the following way. Before and after each glob- 
ular cluster we observed the Galactic A star HIP 2860. Next, 
we convolved a high S/N Mauna Kea night sky spectrum to 
the instrumental resolution and scaled the absorption to match 
strong observed atmospheric absorption lines in HIP 2860. Fi- 
nally, for each globular cluster we interpolated the two scaled 



0.4 - 



0.3 - 



^0.2h 



0.1 - 



10.2 - 



^ 10 

0) 

o 



o 



9.8 - 



9.6 



B189 



j I I I I I I I I I I I I I I I I I I I I I I I [ 



B189 



• Globulor clusters 

• Virgo ellipticols 
■ Como ellipticals 



_L 



_L 



-0.2 -0.1 
[Fe/H] 



0.1 0.2 



Fig. 2. — Comparison of the iron abundances [Fe/H], a- 
enhancement [a/Fe], and ages of the M31 globular clusters to those 
of the elliptical galaxies of vDC. The M31 clusters have similar 
abundance patterns and ages as the elliptical galaxies. Typical er- 
rors in individual measurements are 0.1 -0.2dex, i.e., most of the 
data points are within la of one another. 



absorption spectra obtained from the HIP 2860 observations 
bracketing that cluster This procedure worked well, leaving 
no detectable residuals of night sky absorption lines in the 
wavelength regions of interest. 

One-dimensional spectra were extracted by summing the 
central 13 rows of the two-dimensional spectra, correspond- 
ing to l."8 (7 pc, which is approximately the half-light di- 
ameter of the clusters). The results are independent of the 
extraction aperture; as expected (see § 5) we are not able to 
detect mass segregation in these clusters. The spectra of the 
four most luminous globular clusters were de-blueshifted to 
restframe wavelengths, normalized, and averaged in order to 
increase the S/N ratio and diminish the effects of any remain- 
ing systematic detector or night sky residuals. The averaging 
of these four clusters is justified as their ages and metallicities 
are nearly identical (see Fig. |2). In the following the cluster 
B189 will be discussed separately, although we note here that 
including this cluster in the average spectrum has no impact 
on our analysis or conclusions. 
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Fig. 3. — Averaged Keck spectra of the M31 globular clusters B143, B147, B163, and B 193 (thick grey lines). The left panels show the region 
near the Nal doublet and the right panels show the FeH Wing-Ford band region. The bottom panels zoom in on the IMF sensitive features. 
Colored lines show predictions from stellar population synthesis models with different IMFs. Models with "normal" IMFs match the data very 
well over the entire wavelength range. Dwarf-rich IMFs with steep slopes {x ~ -3) — which provide good fits to massive elliptical galaxies — 
do not fit the observed weak Na I doublet and Wing-Ford band of the globular clusters. 



3. THE NA I DOUBLET AND THE WING-FORD BAND 
IN M31 GLOBULAR CLUSTERS 

Figure |3] shows the averaged spectrum of the M31 globular 
clusters in the wavelength regions near the Na I doublet (left 
panel) and the Wing-Ford band (right panel). Note that the 
Na I doublet is resolved and not blended with TiO, as globular 
clusters have much lower velocity dispersions than the ellip- 
tical galaxies discussed in vDC. Following vDC the spectrum 
was normalized by fitting second-order polynomials in the top 
panels (excluding the regions around the features of interest) 
and first-order polynomials in the bottom panels. The typical 
uncertainty in the averaged spectrum is 0.003 A"', as judged 
from the median scatter among the four globular clusters over 
the wavelength range 8200 A - 8400 A. AH the easily visible 
absorption lines are stellar features, not noise. We smoothed 
the spectra slightly to approximate the resolution of the IRTF 
spectral library (see below). 

Colored lines are stellar population synthesis models. The 
models are the same as those in vDC; they are 10 Gyr old. So- 
lar metallicity models based on empirical stellar spectra from 



the Rayner et al. (2009) IRTF library. Models with slightly 
different ages are virtually identical. The models will be pre- 
sented in detail in C. Conroy et al., in preparation. Differ- 
ent colors indicate different dwarf contributions, ranging from 
bottom-light to bottom-heavy IMFs. The data are well fit by 
models with a Kroupa (2001) or bottom-light IMF: the me- 
dian absolute deviation is 0.005 A"' from 8 150 A - 8700 A 
and 0.007 A"' from 9700 A - 10,200 A. This demonstrates 
that our stellar population synthesis model, which is based on 
Solar metallicity stars in the MiUcy Way, provides an adequate 
description of the spectra of globular clusters in M3 1 . 

We have shown in vDC that elliptical galaxies are best fit by 
a bottom-heavy IMF with x = -3 (the red line in Fig. |3]l. The 
globular clusters are poorly fit by such models. In the Nal 
band (8185 A - 9205 A; see vDC) the median absolute devi- 
ation is 0.025 ± 0.005 for the x = -3 IMF compared to 0.004 
for a Kroupa (2001) IMF In the Wing-Ford band (9910 A - 
9930 A) the median absolute deviation is 0.015 ±0.005 for 
x = -3 and 0.003 for Kroupa. 
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Fig. 4. — Top panels: Direct comparison of the globular cluster spectra to the elliptical galaxy spectra. The globular cluster spectra were 
smoothed to match the average velocity dispersion of the elliptical galaxies. The main sample comprises the luminous clusters B143, B147, 
B163 and B193. The cluster B189 is shown separately; note that the S/N ratio of this spectrum is much lower than that of the averaged 
spectrum of the other clusters. Bottom panels: The black spectrum is the difference between the elliptical galaxies and the main sample of 
globular clusters. The Doppler- broadened Nal doublet and the Wing-Ford band are clearly much stronger in the elliptical galaxies than in the 
globular clusters. The green line is the difference between a stellar population synthesis model with an enhanced dwarf population (x = -3) and 
one with a depressed dwarf population (x = -0.8). The green line is a good fit to the black line, which shows that IMF variations can explain 
the spectral differences between globular clusters and massive elliptical galaxies. 



4. DIRECT COMPARISON BETWEEN GLOBULAR 
CLUSTERS AND ELLIPTICAL GALAXIES 

The fact that stellar population synthesis models indicate a 
normal IMF for globular clusters is significant as those same 
models indicated a bottom-heavy IMF for elliptical galaxies. 
We can also compare the globular clusters to the elliptical 
galaxies in a more direct way. In the top panels of Fig. |4] 
we show the averaged cluster spectrum along with the aver- 
aged elliptical galaxy spectrum. The Coma and Virgo spectra 
of vDC were averaged, and the globular cluster spectrum was 
smoothed to the same velocity width as the elliptical galaxies. 
As expected, the averaged globular cluster spectrum is very 
similar to the averaged elliptical galaxy spectrum, except in 
the regions around Na I and the Wing-Ford band. The glob- 
ular cluster B189 is shown separately. Although the absorp- 
tion lines in this cluster are slightly stronger than in the other 
clusters owing to its higher metallicity, the Nal and Wing- 
Ford band do not reach the depth seen in the ellipticals. This 
demonstrates that our conclusions are insensitive to residual 
metallicity or abundance variations between globular clusters 



and elliptical galaxies. 

The bottom panels of Fig. H] show the difference spectrum, 
obtained by subtracting the averaged spectrum of the globular 
clusters from that of the elliptical galaxies. The two most sig- 
nificant features in the difference spectrum are at ^ 8 190 A 
and at ^ 9920 A. In both cases the absorption in the ellipti- 
cals is much stronger than in the globular clusters. The wave- 
lengths of these features correspond to those of the Doppler- 
broadened Na I doublet and the Wing-Ford band respectively, 
demonstrating empirically that these dwarf-sensitive features 
are significantly stronger in massive elliptical galaxies than in 
M31 globular clusters. Expressed as an index difference, the 
Nal index (as defined in vDC) is 0.027 ±0.007 mag stronger 
and the Wing-Ford index is 0.017 ± 0.006 mag stronger, 
where the uncertainties are determined from the scatter among 
the individual objects entering the two spectra. 

The green line in Fig. |4] shows the difference between two 
of the models shown in Fig.[3j a model with a bottom-heavy 
IMF (x = -3, which was the best fit to the elliptical galax- 
ies in vDC) and a model with a bottom-light IMF (x = -0.8, 
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which provides a good fit to the M/L ratios of globular clus- 
ters [Strader et al. 2011]). The models have identical age (10 
Gyr) and Solar metallicity; the only difference is the number 
of low mass stars. The green line is a good fit to the difference 
spectrum, as might have been expected from the analyses in 
vDC and § 3. The model also reproduces the fact that the Ca 
II A8498, A8542, and A8662 lines are weaker in the globular 
clusters; in the stellar population synthesis model this reflects 
the weak Ca II absorption in dwarfs compared to giants. Al- 
though giants dominate the light, for an IMF as steep as x = -3 
there are sufficient dwarfs to slightly depress the Ca lines. The 
relative weakness of the Ca triplet in elliptical galaxies was 
previously discussed by, e.g., Cenarro et al. (2003). We note, 
however, that in our differential analysis this effect is very 
sensitive to small errors in matching the velocity dispersions 
of the two spectra. 

5. DISCUSSION 

We have shown that massive cluster elliptical galaxies have 
much stronger Na I and FeH absorption than old metal-rich 
globular clusters in M31. This constitutes a critical test of 
the analysis in vDC: equally strong absorption in the globu- 
lar clusters would have implied that the enhancement in el- 
liptical galaxies had incorrectly been interpreted as an IMF 
effect, as the low M/L ratios of globular clusters are incon- 
sistent with dwarf-rich IMFs. Although we cannot exclude 
that some unknown stellar population effect is responsible for 
the strong absorption in elliptical galaxies, it is striking how 
well the difference spectrum in Fig. |4]is fit by an enhanced 
dwarf population. Any explanation not involving dwarf stars 
has to account for the enhanced Na I and FeH absorption in 
ellipticals, the suppressed Ca II, and the similarity of TiO, Fe, 
Mg, H/3 and other elements. We note that known abundance 
anomalies in Galactic globular clusters work in the opposite 
direction (Gratton, Sneden, & Carretta 2004). 

If dwarfs are indeed the cause of the differences in ab- 
sorption then Fig. 4 constitutes direct evidence for significant 
variation in the low mass IMF. Quantifying the IMF in ellip- 
tical galaxies using stellar population synthesis models, we 
can confidently rule out dwarf-suppressed IMFs such as those 
proposed by van Dokkum (2008) and Dave (2008) and also 
"standard" Kroupa (2001) or Chabrier (2003) IMFs which ai-e 
appropriate for the Milky Way disk. This conclusion does not 
depend on the detailed present-day stellar mass function of 
globular clusters (which may have been modified by dynam- 
ical effects), as the expected strengths of the Nal line and 
the Wing-Ford band are almost identical for all "light" IMFs 
(Schiavon et al. 2000; Fig. 3). 

Star formation likely proceeded very differently in massive 
elliptical galaxies than in the disks of spiral galaxies. It is now 
thought that the progenitors of massive ellipticals were very 
compact, with average densities > 20 M0 pc"^ inside the ef- 
fective radius (e.g., van Dokkum et al. 2008; Buitrago et al. 
2008). These densities are similar to giant molecular clouds 
in the Milky Way, but given the ^ kpc scale of these galaxies 
the column densities would have been several orders of mag- 
nitude higher. It will be interesting to measure the physical 
conditions in the star-forming progenitors of these galaxies 
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and compare them to measurements in star-forming regions 
in the Milky Way (see, e.g., Ivison et al. 2010). 

This study can be extended in many ways. The compari- 
son between globular clusters and elliptical galaxies itself is 
model-independent, but we still rely on stellar population syn- 
thesis models to quantify how steep the IMF is in elliptical 
galaxies. The main uncertainty in these models is whether 
short-lived stellar evolution phases with unusual abundance 
patterns are missed, and this can be addressed by augment- 
ing the spectral library. Radial gradients in Na I and FeH can 
provide information on the spatial distribution of the dwarfs, 
which in turn might reflect different formation mechanisms 
for the central of ellipticals and their outskirts (e.g., Oser 
et al. 2010; van Dokkum et al. 2010). Some studies find 
steep gradients in Nal (Boroson & Thompson 1991), which 
might suggest that such effects could be important. It will 
also be interesting to extend this study to elliptical galaxies of 
lower luminosity. There is good evidence that the dynamical 
M/L ratio of early-type systems scales with velocity disper- 
sion (e.g., Cappellari et al. 2006; Forbes et al. 2008; Treu 
et al. 2010). A varying IMF may be responsible for this trend 
(Treu et al. 2010, Dutton et al. 2010), although some stud- 
ies consider dark matter variations a more likely possibility 
(Graves & Faber 2010). As noted by Cappellari et al. (2006) 
the dynamical M/L ratios of low luminosity ellipticals may 
require a normal IMF, as they appear to be lower than the 
M/L ratio of a stellar population with a bottom-heavy IMF 
{M/Lb ~ 23 and M/Lr « 13 for jc = -3). A prehminary anal- 
ysis of a Keck spectrum of the elliptical galaxy NGC 4458 
(cr = 85 kms"'; Cappellari et al. 2006) shows that its dwarf- 
sensitive features are indeed similar to those of the M3 1 glob- 
ular clusters. Finally, it will be interesting to re-examine stud- 
ies of the evolution of scaling relations of early-type galaxies 
(e.g., Renzini 2006). The color and M/L evolution of massive 
galaxies in clusters seem to require a bottom-light IMF with 
X = \ (van Dokkum 2008), but it may be possible to fit 

the data with steeper mass functions when structural evolution 
is incorporated in the modeling (e.g., van der Wei et al. 2008; 
van Dokkum 2008; Hopkins et al. 2009; Holden et al. 2010). 

A varying IMF has important implications for the interpre- 
tation of observations of distant galaxies. The M/L ratios of 
stellar populations with an x = -3 IMF and a Kroupa IMF dif- 
fer by approximately a factor of 3-4. If it is not known how 
the IMF varies with galaxy type, this uncertainty directly cor- 
responds to the systematic uncertainty in the stellar masses 
and star formation rates of galaxies. Even though massive el- 
liptical galaxies make up only a small fraction of the galaxy 
population today, their old ages imply that their progenitors 
become increasingly more dominant at higher redshifts. As 
a result, the star formation history of the Universe, the evolu- 
tion of the galaxy mass function, and many other observations 
would have to be revised. 
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